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^ ' In this work, we introduce two spatio-temporal colored bounded noises, based on the zero- 

dimensional Cai-Lin and Tsallis-Borland noises. 

We then study and characterize the dependence of the defined bounded noises on both a temporal 
CO ' correlation parameter r and on a spatial coupling parameter A. The boundedness of these noises has 

Cn \ some consequences on their equilibrium distributions. Indeed in some cases varying A may induce 

a transition of the distribution of the noise from bimodality to unimodality. 

With the aim to study the role played by bounded noises on nonlinear dynamical systems, we inves- 
fj I tigate the behavior of the real Ginzburg-Landau time-varying model additively perturbed by such 

^ , noises. The observed phase transitions phenomenology is quite different from the one observed when 

a' the perturbations are unbounded. 

. ' In particular, we observed an inverse "order-to-disorder" transition, and a re-entrant transition, 

■y , with dependence on the specific type of bounded noise. 

C^ 
-)— » 

C/3 

^ • Since the hallmark seminal works on Brownian motion by Einstein and Langevin, Gaussian noises (GNs) have been 

i-^ ' one of the main concepts used in non-equilibrium statistical physics, and one of the main tools of its applications, 

^ from engineering to biology. The later, and quite dichotomic, mathematical works by Ito and Stratonovich laid a firm 

O theoretical basis for the mathematical theory of stochastic differential equations, as well as a long-lasting controversy 

, ^, , on which of the two approaches is the best suited in order to better describe stochastic fluctuations in the real world. 

Other hallmarks in stochastic physics have been introduced in the seventies with the birth, in the framework of 

CS| ! the Ilya Prigogine school, of the theory of noise-induced transitions by Horsthemke and Lefever— and in the early 

K*" ■ eighties, in the context of the Rome school, with the introduction of the concept of stochastic resonance-. The study 

^^ ' of stochastic perturbations affecting excitable systems revealed the existence of another stochastic counter-intuitive 

'■ ■ effect: the coherence-resonance^. Finally, in last 20 years nonlinear analysts developed a a rigorous theory of stochastic 

bifurcations both phenomenological and dynamics^. 

Both stochastic resonance and coherence-resonance are example of noise-induced emergence of ordered structures 

^1^ ' in zero- dimensional systems^i^. In other words, noise in a nonlinear system is able to create order. This concept 

Z^ ' stimulated a large amount of research on constructive stochastic effects in spatially extended systems. A first important 

result of these investigations has been the emergence of noise-induced phase transitions " disorder-to-order" in spatially 

extended systems. In those cases a key bifurcation parameter is the strength of the spatial coupling (e.g. in the case of 

Laplace coupling, the diffusion coefficient). It is remarkable that this phenomenology was also observed in the case of 

stochastic perturbations of deterministic physical systems unable to form patterns, i.e. when moving from a spatially 

^ homogeneous to a spatially coupled deterministic system there is no emergence of Turing patterns. The addition of 

^_' noise to such systems may cause the onset of "disorder-to-order" phase transitions. An example of such a model is 

the Ginzburg-Landau (GL) real equation^"—. 

The majority of works in the field of noise-induced phase transitions were based on white noises, in which is absent any 
spatio-temporal correlation or structure. This kind of fluctuations is appropriate when modeling internal "hidden" 
degrees of freedom, of microscopic dynamical nature. On the contrary, fluctuations originated externally to the system 
in study may exhibit both temporal and spatial structure -iS. Zero-dimensional systems perturbed by colored Ornstein- 
Uhlcnbeck noise (OU) exhibit correlation-dependent properties that arc missing in case of null autocorrelation. Among 
them: the emergence of stochastic resonance also for linear systems, and re-entrance phenomena, i.e. transitions from 
monostability to bistability and back to monostabilitj*i2rJ^. Spatially extended systems exhibit even more striking 
effects when they are perturbed by spatially white but temporally colored noises, due to a complex interplay between 
noise intensity, spatial coupling (of the perturbed system) and autocorrelation time^°. 
Garcia- Ojalvo, Sancho and Ramirez-Piscina introduced in^'^ the spatial version of the Ornstein-Uhlenbeck noise. 



characterized by both a temporal scale r and by a spatial scale A, and studied analytically the spatio-temporal 
correlation of this kind of noise io^. Later they showed that the interplay of this noise with suitable nonlinear 
systems can induce transitions often unobserved when the applied noise is white. In particular, when analyzing 
the additive spatio-temporally colored perturbation of the Ginzburg-Landau field, they showed the existence of a 
nonequilibrium phase transition controlled by both the correlation time and the correlation lengthSii^. Similarly to 
the deterministic bifurcations ruled by two parameters, the transitions induced by the Garcia-Ojalvo, Sancho and 
Ramirez-Piscina spatio-temporal noise (GSR) show a more complex picture than the transitions induced by a spatio- 
temporal noise that is spatially uncorrelated^^. In the case of a single spatial dimension, a different extension of the 
Ornstein-Uhlenbeck noise was proposed, and analytically studied ini^. 

The above-summarized body of research is essentially based on the use of GNs, whose background is the Central 
Limit Theorem, and which is the best approximation of reality in many cases. However, an increasing number of 
experimental data motivated theoretical studies stressing that many real-life stochastic processes does not follow white 
or colored Gaussian laws, but other densities (such as fat-tail power-laws). More recently, a vast body of research 
focused on another important class of non-Gaussian stochastic processes: the bounded noises. The studies on bounded 
noises, apart some sporadic exception, were mainly confined to the Dichotomous Markov Noise (DMN) and to its 
applications-. In the last twenty years, together with a renewal of theoretical interest for DMN^, other classes of 
bounded noises were defined and intensively studied in statistical physics and in engineering, and - to a lesser degree 
- in mathematics and quantitative biology. 

The rise of scientific interest on bounded noises is motivated by the fact that in many applications both GNs and fat- 
tailed non-Gaussian stochastic processes are an inadequate mathematical model of the physical world because of their 
infinite domain. This should preclude their use to model stochastic fluctuations affecting parameters of dynamical 
systems, which must be bounded by physical constraints. Moreover, in many relevant cases, especially in biology, 
some parameters must also be strictly positive. As a consequence, not taking into account the bounded nature of 
stochastic fluctuations may lead to unrealistic inferences. For instance, when the onset of noise-induced transitions 
depends on exceeding a threshold by the variance of a GN, this often means making negative or excessively large a 
parameter. To give an example taken from medicine, a GN-based modeling of the unavoidable fluctuations affecting 
the pharmacokinetics of an antitumor drug delivered by means of continuous infusion leads to a paradox. Indeed, the 
probability that the drug increases the number of tumor cells may become nonzero, which is absurd^^. Thus in order 
to avoid these problems, the stochastic models should in these cases be built on bounded noises. 
The deepening and development of theoretical studies on bounded noises led to the attention of a vast readership on 
new phenomena, such as the dependence of the transitions on the specific model of noise that has been adopte d^^'^^ . 
This means that, in absence of experimental data on the density and spectrum of the stochastic fluctuations for the 
problem in study, a scientific work should compare multiple kinds of possible stochastic perturbations. Moreover, 
currently the bounded noise approach also implies that the possibility of obtaining analytical results is remarkably 
reduced. For example, the study of stochastic bounded perturbations acting on a simple scalar deterministic model 
leads to at least two stochastic differential equations. Indeed, one or more additional equations must be devoted to 
the modeling of the bounded stochastic processes. 

Finally, we remind that in order to generate a temporal bounded noise, two basic recipes have been adopted so far. 
The first consists in applying a bounded function to a random walk^^, whereas the second one consists in generating 
the noise by means of an appropriate stochastic differential equatio n^^i^^ . 

Our aim here is twofold. First, we want to define two simple families of spatio-temporally bounded noises, which 
extend two kinds of temporal noises frequently employed in literature, the Tsallis-Borland noiseSS, and the Cai-Lin 
noiseSi. 

Second, we want to illustrate the possible effects of external bounded stochastic forces (i.e. of additive bounded noises) 
on nonlinear spatio-temporal models. As a case study we have chosen the scalar Ginzburg-Landau equation, one of 
the best-studied amplitude equation representing "universal" nonlinear mechanisms. 

Phase transitions induced in GL model by additive and multiplicative unbounded noises were extensively studied for 
this model^i^ i^^i^^" — . GL equation driven by temporally varying Markov dichotomous noise have been considered 
in3^. Moreover, recently it has been investigated the case of GL dynamics under external additive "quenched noise" 
(i.e. a random perturbation that is very slow in comparison to the unperturbed model dynamics) has been studied 
in^, where the possibility of multistability was analytically studied. 

Simulations with the GSR unbounded noise will be performed, in order to dissect the effects that are more closely 
related to the boundedness of the noises, and to their type. 

Among the tools used for the analysis of the numerical simulations, we shall also use the Moran mde:i^^—, a classical 
tool of spatial statistics of lattice data, which has been intensively applied in geographical statistics^ and in image 
processing^!, which - however - has never been used in statistical physics, at the best of our knowledge. 



II. BACKGROUND ON TEMPORAL BOUNDED NOISE 

Here we provide some background material related to the generation of temporally correlated bounded noises. 

A. The Tsallis-Borland bounded noise 

The first noise we take in consideration is the Tsallis-Borland bounded nois o^"'^^ , whose dynamics is described by 
the following equation: 

A Tsallis q-statistics22ii2^ is the stationary distribution of this noise: 

PTB{O^A{B^-e)t'\ (2) 

whose mean is null, and whose standard deviation is given by 

CrTs(g) = B J- 

2r(| + T^ 

The parameter q G [—00, 1] links the true autocorrelation time Tc of ^(i) to the parameter r as follows^S: 

5-3g 

T « Tr , 

2 

while the coefficient D is in turn linked to the bound B and to the Tsallis parameter q by the following relationship: 



'2D^^I2 1-q 



TC 5-3q 
In the Supplementary Materials (SMs) we define a generalization of the Tsallis-Borland noise. 

B. The Cai-Lin bounded noise 

j- ^21^33,34 ^ the following family of bounded noises was introduced: 

eit)^^-m+9{0vit), (3) 

Tc 

with 5(|-B|) = 0. The bounded nature of the noise described in |3] easily follows from the fact that at ^ = +B it is 

^' < 0, whereas at ^ = -B it is £,' > 0. 

Note that if g(C) is symmetric then the process ^(t) has zero mean, and the same autocorrelation of the OU process2ii2^, 

i.e. Tc denotes the real autocorrelation time of the process £^{t). 

Observe that, as shown in^-'^''^^, a pre-assigned generic stationary P{£,) can be obtained by the eq. ([3|). Indeed, by the 

Fokker-Planck equation associated to eq. ([3|) it is straightforward to verify that the following function g{£,): 



"^'^f^/l"^'"'""- 



has to be chosen in order to recover P(^). For example, by choosing the following symmetric functior 



21,33. 
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one gets the following stationary density: 



PcL{0-A{B'-e)% (5) 



parametrized by S > —1, whose standard deviation is given by: 

r((5 
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For (5 > the distribution is unimodal and centered in 0, while for —1 < (5 < it is bimodal, having a "horns" -like 
distribution with two vertical asymptotes at ^ — ^ ±B. 

III. SPATIO-TEMPORAL COLORED NOISE: BACKGROUND AND NEW DEFINITIONS 

Let us consider the well-known Ornstein-Uhlenbeck equation: 

e{t)^ — m + ^^vit), (6) 

where ri(t) is a Gaussian white noise of unitary intensity: (r7(i)r7(i — ti)) = d{t — ti), and whose solution is a Gaussian 
colored stochastic process with autocorrelation: 

{ri{t)ri{t ~ h)) ^ A exp ' '^"*'' 

Equation ^ has been generalized in2^ as follows: 

1 V2D 

dtax,t)^ ax,t) + ry(a:,i), (7) 

Tc Tc 

where 'r]{x,t) is a noise white in space and time. As a consequence £,{x,t) is a stochastic process white in space but 
colored in the temporal dimension. Finally, a noise also colored in space was proposed in22, by adding to eq. ([7]) the 
most known and simple spatial coupling, the Laplace operator, thus yielding the following partial differential Langevin 
equation 

2^2 ^ \/2D 

dt^x, t) V2^(a:, t) - -^(x, t) + 7^{x, t), (8) 

Tc Tc Tc 

where A > is the spatial correlation lengthi^ of ^{x,t). As usual in non-equilibrium statistical physics, in line 
witbi^iSij we shall investigate the lattice version of eq. ([8]), i.e.: 

\2 1 \/9D 

^pit) = TT^lUt) -Ut) + ^^Vpit), (9) 

^Tc Tc Tc 

where p = h x{i,j) is a point on a M * M lattice with steps equal to h. The symbol V| denotes the discrete version 
of the Laplace's operator: 

^lUt)^h-' E (<^'-<^p)' (10) 

iGNe{p) 

where Ne{p) is the set of the neighbors of the lattice point p. 

In order to define spatio-temporal bounded noises based on the Tsallis-Borland and on the Cai-Lin noises, we shall 

use an approach analogous to the one employed i n^^i^^ to extend the OU process. 

Thus the Tsallis noise can be generalized to spatially extended systems by means of the equation: 

dta^,t) = — V2^(x,i) - -—^-L— + ^^rj{x,t), (11) 

2tc t1~ {^{t)/By T 

while the spatio-temporal Cai-Lin noise can be defined as follows: 



o\2 -1 / or) 

dta^^t) = —v'a^,t) ~ -e(.T,t) + J y w^rj{x,t). (12) 

Tc Tc \j Tc[l+d) 

We shall study the above two equations in their lattice-based discretizations. 



IV. STATISTICAL FEATURES OF SPATIO-TEMPORAL NOISE 



We mainly characterized the global behavior of the bounded nosies defined by the lattice versions of eqs. (jlip and 
(I12p by means of the well known order parameter— 
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and of its relative fluctuations cr^ (S) , defined as^. 



a^(S)^ 
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Concerning, instead, the characterization of the spatial-autocorrelation-related properties of the spatio-temporal noises 
treated in section IIIIl and in the study of phase transitions induced by them, we employed a classical tool of spatial 
statistics: the Moran's inde x^^i'^° . Given a field ai on a lattice, let us define the following field Xi on the same lattice: 



^ ' ieNe{i) 



ai = +ai + h^V\ai, 



the Moran's statistic In is defined as follows: 



/„-PC(a,x(a)) 



(13) 



(14) 



where PC{u,v) is the Pearson correlation between two vectors u and v. Moreover, in order to compare the degree of 
mutual correlation of two different lattice fields a and 7, some years ago we defined a bivariate Moran inde»^: 



Ir. 



PC{x{a.),xh*)), 



(15) 



where the * in the pedix denotes the normalization of the lattice field, e.g.: a* = (a — {a})/ y^Var{a). 
We first investigated the dependence of the bounded spatio-temporal noises defined by eq. (|TT|) - (fT2|) on the parameters 
Tc (measuring the temporal autocorrelation) and A (measuring the spatial coupling). For both Tsallis-Borland and 
Cai-Lin noises, in our simulations we obtained that the fluctuation a^ is independent of both tc and A (see fig. [TJa 
and l.a of SM). On the contrary, as expected, the Moran index strongly depends on A, and it is independent on tc 
(see fig. [TJb and l.b of SM). We also simulated the GSR unbounded noise by empirically setting (in order to get a 
noise strength somewhat comparable to the one employed in the bounded cases) B = 2(7„, where au = y/2D is the 
strength of the GSR noise, i.e. the standard deviation of the generative white noise in equation ©. We obtained that 
the fluctuation parameter a^ in this case is quite large and it is -differently from the bounded noise case - a decreasing 
function of t (see fig. 2 of SM). Instead, the Moran index behaves as in the previous case. 
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FIG. 1: Tsallis-Borland bounded noise : characterization of spatio-temporal fluctuations. Panel (a) plot of cr(H); panel (b): 
Moran index. Parameters: 40 x 40 lattice, B = 1 and q = —1. 



The global behavior of the spatio-temporal noise, analyzed by the macroscopic observable S, has its counterpart in the 
equilibrium heuristic probability density of the lattice variables ^p. First, we observed that in both types of bounded 
noises, the distribution of ^p is independent of the temporal correlation parameter (see figures 3 and 4.b of SM). On 
the contrary, the spatial coupling parameter A deeply affects the distribution of S,p, and in a noise- type dependent 
manner. Indeed: 

• for Cai noise with 6 < 0, the increase of A induces a transition from a bimodal density to a unimodal density 
(see panel B of figure [5]) ; 

• for both Cai noise with S > and for Tsallis noise, the increase of the spatial coupling induces a decrease of the 
variance of the noise (see fig. [2] and figure 4. a of SM). 
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FIG. 2: Influence of cha racteristic correlation length A on Cai bounde d no ise equilibrium distribution Peq{^p) of a 40 x 40 
lattice system. Panel (a) parameters B = 1, t = 2 and S — +0.5. Panel (b) same noise with t — 1 and 5 — —0.5. 



BOUNDED-NOISE-INDUCED PHASE TRANSITIONS IN TIME-DEPENDENT 
GINZBURG-LANDAU MODEL 



In this section we compare the classical Gaussian approach to noise-induced transitions, to the new approach 
based on bounded noises. In particular, we focused on a well-defined case study in Gaussian noise-induced phase 
transitions?^^: the real Ginzburg-Landau stochastically perturbed equation. We remark here that, although we shall 
show some effects of interest, here our aim is simply to illustrate the differences in the behavior of this system in 
response to bounded and unbounded noises. 
jjj8ji5 Garcia-Ojalvo et al. studied on a two dimensional lattice the real Ginzburg-Landau model: 



1 



Mp = 7;{i'p-i'l + vi Vp) + ^(x, t) 



(16) 



with external noise ^ given by equation ([9]). The worka^^i^ showed that botli spatial and temporal characteristic 
lengths of the noise ^(a;, t) promote the order phase in the face of the Gaussian spatitoemporal noise ry, i.e. they shift 
the order/disorder transition point on the right. 

In the following we present some results from our extensive numerical study of GL equation by modeling the additive 
noise appearing in eq. (IT51) with Gaussian, Cai-Lin and Tsallis-Borland bounded noises. Note that in order to properly 
compare the different behaviors of the system with bounded and unbounded noise, we set:aunbound = (Abound- 
Similarly to the study of spatio-temporal noise described in section IIVI we characterized the phase of GL model by 
means of the global magnetization M =< | '^u ,) V'i.jl > \-^^! of its relative fluctuations (j\[, as well as by using both 
the Moran index of the field tpp (/^) and the bivariate Moran index between the noise field ^p and the GL one tpp 

In our numerical study, the initial condition of GL system was the ordered phase, i.e. '0p(O) = 1, thus we measured 
the robustness of order in the face of noise. 



In absence of spatial coupling in the noise (A = 0), both in the Cai and in the Tsallis case the behavior of the 
magnetization M when varying the autocorrelation time of the noise suggests that bounded noise seems to promote 
the disorder phase for GL system. Indeed, for bounded noises the magnetization is a decreasing function of Tc, see figure 
[31a (and figure 5 of SM). Thus temporal correlation has the opposite effect with respect to the unbounded GSR noise, 
which enhances the order phase^^. It is possible to roughly explain this effect by considering the different behavior 
of the equilibrium distributions Peq{^p) in the bounded and in the unbounded noise. Indeed, in the unbounded noise 
the standard deviation of f scales with 1/tc, and the transition from disordered to ordered phase with Tc could be 
considered as an effect of noise amplitude reduction. On the contrary, both in Cai-Lin and Tsallis-Borland noises 
the equilibrium distributions do not change their standard deviation, as one can deduct by PeqiC) formulas ([2]) and 
([S]), which are independent on r^. In addition, the numerical study of the bounded noise equilibrium distribution for 
A ^ 0, reported in the SMs (figure 3 and 4), suggests that the influence of bounded noise in order/disorder transition 
would be the same also in presence of a fixed spatial correlation. Thus here the overall effect of Tc consists only in 
increasing the spatial autocorrelation of the noise. As a consequence, the field "0 is driven by an even more quenched 
noise, with a constant broad distribution, which enhances the disorder phase. 

For Cai-Lin noise with (5 > 0, the qualitative behavior of ctm and of the Moran index /* are similar (see fig. [Hb 
and[3]c): they both show a maximum at transition point, for Tc ~ 0.5. In the case of Tsallis-Borland noise, M, gm 
and the Moran index are similar to the ones obtained in case of Cai-Lin noise with 5 > 0, suggesting that the GL 
phenomenology with temporally correlated bounded noise is similar in case of unimodal noise distribution. 
In the case of Cai-Lin noise with bimodal distribution (5 < 0) the behavior of GL equation around the transition region 
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FIG. 3: Absence of spatial coupling in the noises (A = 0). Effects of temporal autocorrelation Tc on GL model perturbed by 
additive spatially uncorrelated noises. Noise used: unbounded (GSR) and Cai with 5 = ±0.5. Top panels: global magnetization 
M |(a)| with its relative fluctuations ctm (b) Bottom panels: Moran index of field ^ |(b)| and bivariate Moran index Jj^ between 
the field {■(/;} and its corresponding noise one {^} (d) Parameters: 40 x 40 lattice, A = 0. Cai noises, with B = 2.4, are compared 
with a corresponding unbounded noise with same standard deviation, i.e. y/2D = BjA = 0.6 for 5 = 0.5 and y/2D = B/2 = 1.2 
for 5 = -0.5. 
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FIG. 4: Effects of Cai parameter S on the order/disorder/order A transition. From the observation of M one may gather that 
the first ordered phase is enhanced by the steepness of the bimodal distribution of the noise, which tunes both the level of 
magnetization and the critical value of A. Parameters; 40 x 40 lattice, B — 2.6 and r — 0.3. 



is different (see again fig. [3]b and[3lc): ctm is a decreasing function of Tc, whereas the Moran index steadily increases 
with the temporal autocorrelation. Actually, our numerical results support the idea that this order/disorder transition 
is smoother than a second order one. In a coarse-grained description one could say that the "quasi" dichotomous 
feature of the equilibrium distribution allows the noise to jump back and forth between two symmetric values. Thus, 
in each GL field point the temporal mean value of noise is very close to zero. On the contrary, if the distribution is 
unimodal, the dispersion term slowly approach the noise to the central value 0, so that the permanence on the same 
sign is longer and the temporal mean of noise in a lattice point is greater. Bimodal noise is more prone to ordered 
state, resulting in an increasing order in the ordered phase and in a smoother transition. 

If one equates the variance of the bounded noise to the one of the Gaussian noise, the bound B results to be larger than 
the standard deviation of the correspondent Gaussian noise. In spite of this, for low Tc our two classes of bounded 
noises order the system better than the unbounded noise, enforcing the idea that rare value of unbounded noise, 
placed at the tail of distribution, could have a non-trivial role in phase ordering phenomena. 

Finally, by computing the bivariate Moran indexes I^^ between the discretized field {V'p} and the noise {£,p}, we 
obtained that (see fig. [Hd): i) unbounded noise is not correlated or, for large values of Tc, poorly correlated with the 
field; ii) for bounded noises, the index I^^ is not null, and it increases with Tc, supporting the idea that the net effect 
of temporal correlations consists in freezing the noise, giving it quenched features. In the unbounded noise this effect 
is hidden by the reduction of noise amplitude, which drives the disorder/ordered transition. 

Summarizing, in absence of spatial coupling in the bounded noise, the transition is mainly ruled by the temporal 
autocorrelation of the noise, i.e. by its level of "quenching". This behavior is at variance with the one corresponding 
to unbounded perturbations, where the transition is driven by the amplitude of noise. 

The behavior of the system is deeply affected by the presence of spatial coupling in the noise. Note that we assessed 
the infiuence of the coupling parameter A in correspondence of a relatively low noise temporal autocorrelation r^ = 0.3, 
so that the noise cannot be considered quenched-like. Indeed, the characteristic time of evolution for small ip in the 
deterministic case £,p{t) = is 2. If the noise is of Cai-Lin type (see figure |4] and [5]) , a deep difference in the behavior 
of the magnetization in function of A emerges for different values of 6, as illustrated in fig. 21 With S > 0, similarly to 
the response to unbounded noise, there is a transition towards order, very likely due to the reduction of the variance of 
the distribution of the ^p's. On the contrary, in the case S = —0.5 there is a re-entrant transition order/disorder/order. 
This is very likely a consequence, although non-trivial, of the stochastic bifurcation in the distribution of the ^p's, 
which passes from bimodality to unimodality close to the point of minimum of M. For S — —0.75 the magnetization 
has two extremal points. 

For d ~ —0.5 the spatial autocorrelation of V', measured by /^, rapidly increases in that critical region; for 8 — 0.5 
the Moran index initially decreases with A (which is quite counter- intuitive) and then increases (fig. [5]c). Finally, the 
bivariate Moran index has a pattern that is on the whole increasing for 5 — -t-0.5, whereas for 8 = —0.5 it exhibits a 
minimum (fig. Od). 

In the case where the additive perturbation term in the GL equation is a Tsallis-Borland noise, one may again observe 
a re-entrant transition, although somewhat different from the corresponding transition found for Cai-Lin noise with 
same noise equilibrium distribution. This can be deduced by comparing figures [3] and [BJ where magnetization of GL 
field is shown for three values of the Tsallis parameter q. 

Moreover, figures |4] and |6] suggest that the second ordered phase is enhanced by reducing the standard deviation of 
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FIG. 5: Effects of the spatial correlation parameter A in case of unbounded GSR and Cai-Lin noises with S = ±0.5 in the 
additively perturbed GL model. Top panels: global magnetization M (left) with its temporal variance aM (right). Bottom 
panels: Moran index of field ^p (left) and bivariate Moran index 7{^ between the field ^/j and its corresponding noise ^ (right). 
Cai-Lin noises, with B = 2.6, are compared with a corresponding unbounded noise with same standard deviation, i.e. \/2D — 
B/4 = 0.65 for S = 0.5 and V2D = B/2 = 1.3 for S = -0.5. Parameters: 40 x 40 lattice, r = 0.3. Unbounded noise has 
V2D — 1.3, while Cai noises has B = 2.6. 

the noise distribution, i.e. increasing 6 and q. This consideration supports the idea that the second ordered phase 
is obtained through a mechanism of noise amplitude reduction, similarly to the unbounded case. On the contrary, 
the first ordered phase emerges for reasons more intrinsically related to the interplay between the boundedncss of the 
distribution and the spatio-temporal features of noise field. However, the mechanisms that allow for this phase are 
different in Cai-Lin and Tsallis-Borland noise, because the physics that governs their dynamics is different. 
In all the above-examined cases, the comparisons of the different statistics made it possible to qualitatively discriminate 
between the two ordered phases observed in the re-entrant transitions. Indeed, the first one, which appears at small 
A, is characterized by low values of I^ (blue curve in fig. [5jc, for A < 0.5), and in it order is enhanced by the low 
temporal correlation of the noise. On the contrary, the second ordered phase, appearing at larger A, is characterized 
by large values of /^ (fig. Elc, for A > 1.5), and in it the order is caused by both the lower amplitude of the noise 
and by the strong spatial correlation of the noise. Finally, in figure |S] we show the effect of the amplitude B on the 
distribution of the GL field. Increasing B, for non-small values of A one moves from unimodal density mainly located 
close to '0 = 1 to a bimodal density. 



VI. CONCLUDING REMARKS 



In the first part of this work, we defined two classes of spatio-temporal colored bounded noises, which can be 
directly derived by two zero-dimensional bounded noises: the Tsallis noise and the Cai noise. We have analyzed the 
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FIG. 6: EfTects of Tsallis parameter q on the order/disorder/order A transition. From the observation of M one may gather 
that the region occupied by the second ordered phase increases with q, i.e. inversely with the variance of the noise distribution, 
which tunes both the level of magnetization and the critical value of A. Parameters: 40 x 40 lattice, B — 2.6 and r — 0.3. 
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FIG. 7; Effects of the spatial correlation parameter A in case of unbounded and Tsallis-Borland noises with g = — 1 in the 
additively perturbed GL model. Panel (a) global magnetization M with relative fluctuations aM- Panel (b) Moran index of 
field ^p and bivariate Moran index I^^ between the field {^} and its corresponding noise one {^j. Tsallis noise, with B — 2.6, 
is compared with an unbounded case with same standard deviation, i.e. \/2D — B/4 = 0.65. With bounded noise, GL system 
shows an double order/disorder/order transition. 



role of the spatial coupling parameter A and of the temporal correlation parameter r on the distribution of the noise 

by employing suitable statistical observables. Differently from unbounded noise case, the equilibrium distribution of 

the studied bounded noises does not depend on r, while in some cases the increase of A induces transitions from 

bimodality to unimodality in the distribution. These features could be relevant when bounded noises are applied to 

dynamical systems, in particular in the presence of noise-induced phase transitions. 

With this in mind, in the second part we employed the above-mentioned two kinds of bounded noises to investigate 

the phase transitions of the Ginzburg-Landau model under additive stochastic perturbations. 

Our simulations showed a phenomenology quite different from the one induced by colored unbounded noises. 

To start, in presence of spatially uncoupled bounded noises, the increase of the temporal correlations enhances the 

quenching of the noise, eventually producing an "order-to-disorder" transition in the GL model. Note that when the 

perturbation is unbounded an opposite transition is observed. 

Furthermore, spatial coupling induces contrasting effects on the spatio-temporal fluctuations of the noise, resulting for 

some kind of noises in a re-entrant transition (order/disorder/order) in the response of the GL system. This specific 

case of noise type-dependence in spatio-temporal dynamical systems is novel, and is in line with previous observations 

in zero-dimensional systems"'^^'^-. 

Finally, we stress that in defining the spatial extensions of the Cai-Lin and of the Tsallis-Borland temporal noises we 

adopted as spatial coupling the classical Laplace operator, in line witfti^. However, more sophisticated alternatives 
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FIG. 8: Equilibrium distribution Peq{'4^p) of GL system on a 40 x 40 lattice, for Cai noise with different sizes of bound. Noise 
parameters are S = —0.5 and r = 0.3. 

might be employed, to take into account long-range interplay a^'^^ provided that the constraint of boundedness of the 
noise can be preserved. 
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CN ■ I. A GENERALIZATION OF THE TSALLIS-BORLAND BOUNDED NOISE 

T — I 

o . 

fvq , The Tsallis-Borland bounded noise noise is here generalized as follows: 



i3: ^'{t)^-f{o + ^^vit), (1) 



en 
(N 



(N 



X 



1 „,.s V2D 



T T 



where f{+B) = — cx) f{—B) = +oo. Our extension of the Tsallis-Borland noise allows to generate a noise with a 
pre-assigned stationary density P, . Indeed, from the Fokker-Planck equation associated to eq. ([l} , one immediately 
obtains that the following /(^) must be chosen to be : 



^'. DP' (e) 

O. m = ^^y (2) 

^ For example, in order to obtain the density P{£,) = (l/2)Cos(^)+, one must consider: 
-(— > 

C3 , D 

^ . /(O - --taniO. 

B 

^ I Note that in Tsallis noise the unboundedness of /(^) at ^ = |i?| may cause problematic stiffness during numerical 

Q . simulations. For example in the Euler schema 

O 



II. USEFUL TRANSFORMATIONS 



1 \/2D 

C(i + At) = ^{t) + - fiOAt + Gt\/Ai 

r r 



(N 

Q if £,{t) is close to \B\ then it may be \^{t + At)| > B. 

f — . A possible remedy, suggested for the Tsallis-Borland noise in ref. 20 of the main text (H. S. Wio and R. Toral. Physica 

PsJ D, 193, 161-168, 2004) is to re-extract the gaussian random variable Gt- 

l/~) Another trick, which we introduce here, consists in applying the following nonlinear transformation to 

g : ^ = B szn(<^), (3) 



which yields the following Langevin equation for 0(i): 

c/,' = a(0) + 6(</))77(i), (4) 

where: 

tB cos[(p) 

/ ,x ^ D sin{4>) f{B sinjcp)) 
"^'^^ B^T^cos^q^) Bcos{cj)) ■ ^^ 

Thus, in the case of temporal Tsallis noise (equation (1) in the text) we have: 

1 + q sinjcj)) 

aw) = ^ J7-T (7) 

2t cos'^[(pj 



In the spatiotemporal case, similarly to the purely temporal case it is useful to define the transformation: 

ip = B sin{4>p) 
leading to the N equations: 

0;(i)=ap(</))+6p(./))77p(t). (9) 

where: 

bpW^^ T-T (10) 

, ,, D sin((h„) f(B sin(6r,)) A^ 1,9 v^ , ■ ,^ ^ ■ , , ^^ /,,n 

Namely, in the case of spatiotemporal Tsallis noise (equation (11) in the text) we have: 

fA\ (1 + g) sin{(j)p) A2 1 ^_2 V^ f ■ fA \ ■ fA w no^ 

"^^^^^~ 2r co.3(0,) +^^^;(^^ ^ L^^^ (--(^") - -n(0,)) (12) 



2 M-g ] 1 



bpi^)-\l-\^^.]zZZTT-^- (13) 



Similarly to the Tsallis' noise, also in simulations of the Cai noise one may, although less frequently, encounter 
numerical stiffness. By applying the transformation: 

^ = B sin{(l)) (14) 

one gets a stochastic equation of the form of eq. ([4]), where: 

g{B sinjct))) 
B cos[(p) 

In the case of temporal Cai-Lin noise (equation (3) in the text) we have: 

a(0) = -^(l^*-(^) (17) 



K(t>) = 1 / /-^ _^ g. Signum(Cos{4>)) (18) 



When simulating the spatiotemproal Cai-Lin noise, applying the transformation ([T4| yields the following N stochastic 
equations: 

(/>p(i) = ap{(j)) + bp{(j))r]p{t), 

where: 

g(5 sm(0p)) 



3 

/ I 1 \ A^ 1 

n£Ne{p) 

Namely, in case of the spatiotemporal Cai-Lin noise (equation (12) in the text) we have: 

\ -L 25 A^ 1 

«('/') ^ - ^. -I ^A-i ^Q"('?^p)+9 7T^^^^ Yl {sin{(f)n) - sin{(f>p)) (21) 



^p(0) ^ \/ T- (i I g-. Signum{Cos{(l)p)). (22) 



III. NUMERICAL ALGORITHMS 



Both Cai and Tsalhs spatiotemporal noise generative equations, once reduced to the form ([3]), by mean of the 
trasformation ^, can be considered has a Langevine equation in (j) with multiplicative noise. We have integrated this 
noise numerically with the derivative free Milstein method for Ito rapresentation, i.e.: 



Undt + bnVdtrj, 

dt r 



t>n+l = (t>n + dndt + bnV dtrjn H 



{4>n)- bn [vl 



-I 



Once noise values ^ = Bsm{(j){t)) have been obtained, GL equation is integrated with a standard Runge-Kutta 
method: 



dtiPp = i^ (vJ") + ip 
^p = Mi) + ^ (f [^) + ^p) 
^p = Mt) + f (^ (^) + ^ (V')) + dt^P- 



We have used a time step dt = 0.005, for a total time of T = 250, taking the mean values of the observables in the 
last interval of lenght 25, whit a number of different runs between 10 and 20. 

IV. SUPPLEMENTARY FIGURES 

V. MORAN'S SCATTERPLOTS 

Apart the Moran Index, in order to assess the correlations between a lattice field a and the field x constituted by the 
mean of the values of a in the neighbouring points at each point p of the lattice, we also employed the scatterplot of 
the couples {ap,x{cep)), called Moran scatterplot. This is an useful tool of descriptive spatial statistics. For example, 
in fig. [niwe show the effects, evaluated by means of the Moran's scatterplots, of the temporal autocorrelation Tc on a 
discrctized GL model in a 40 x 40 lattice perturbed by Tsalhs noise. 





(a) 



(b) 



FIG. 1: Characterization of spatiotemporal fluctuations for bounded Cai noise: plots of a (H) (a) and Moran index (b) 
Parameters: 40 x 40 lattice, B = 1 and 5 = —0.5. 





(a) 



(b) 



FIG. 2: Characterization of spatiotemporal fluctuation for unbounded spatiotemporal noise. Plots of a (H) (a) and Moran 



index (b) Parameters: 40 x 40 lattice, V2D = 0.5. 
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FIG. 3: Influence of characterisctic correlation length r on Cai bounded noise equilibrium distribution Peq{S,p) of a 40 x 40 
lattice system. Panel (a) Cai noise with _B = 1, A = 1 and 5 = +0.5. Panel (b) same noise with A = 0.5 and S = —0.5. 
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FIG. 4: Influence of r and A on TsaUis bounded noise equilibrium distribution Peq{S,p) of a 40 x 40 lattice system. Parameters 



are B = I, q = —1, r = 2 for panel (a) and \ = 2 for panel (b) 
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FIG. 5; Effects of temporal autocorrelation r on a discretized GL model in a 40 x 40 lattice perturbed by TsaUis noise, with 
q = —1 and B = 2.4 and A = 0. GL phases are characterized by the global magnetization A/, the relative fluctuations aM, the 
Moran index of field tp and the bivariate Moran index I^^ between the field {^p} and its corresponding noise one {^}. With 
bounded noise, GL system shows an inverted order/disorder transition with respect to the unbounded case. 





FIG. 6: Effect of spatial correlation strength A on Moran scatterplots of the GL model with additive noise 
scatterplots show the effect of the spatial coupling A on the GL model 
transition induced by Cai noise {5 — 



Panel (a) 



The Moran 
scatterplot of the order/disorder/order 
-0.5, B = 2.6); panel[(b)j TsaUis noise (g ='^, B = 2.6). In aU panels r = 0.3. Spatial 
coupling in both cases increases the slope of scatterplot, and at the same time has a tendency to make the cloud of the scatter 
plot more compact, lowering the dispersion of short range correlation around the mean value. 



